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Abstract 

We discuss phenomenology of the vacuum condensate ((A'} l ) 2 nin ) in pure gauge the- 
ories, where A" is the gauge potential. Both Abelian and non-Abelian cases are 
considered. In case of the compact U(l) the non-perturbative part of the condensate 
((^u)min) * s saturated by monopoles. In the non-Abelian case, a two-component 
picture for the condensate is presented according to which finite values of order 
A.qq D are associated both with large and short distances. We obtain a lower bound 
on the by considering its change at the phase transition. Numerically, it 

produces an estimate similar to other measurements. Possible physical manifesta- 
tions of the condensate are discussed. 

1 Introduction 

Vacuum condensates are commonly used to parameterize non-perturbative effects in the- 
ories with strong couplings. Focusing on QCD, the best known vacuum condensate seems 
to be the quark condensate (qq) . Perturbatively, the quark condensate is proportional to 
the bare quark mass, m q and vanishes in the chiral limit m q = 0. Thus, a non-vanishing 
quark condensate signifies spontaneous breaking of the chiral symmetry. 

There are also condensates which do not vanish perturbatively. A well known example 
is the gluon condensate, (ot s {G a ^ v ) 2 ) where is the gluon field strength tensor [[]J. 
Perturbatively, (a s (G^,) 2 ) ~ K\j V where A uv is the ultraviolet cut off. Thus, care should 
be exercised to separate the non-perturbative contribution from the trivial perturbative 
part. There are various ways of subtracting the perturbative contribution. Historically, 
the QCD sum rules were the first example of such a subtraction. Here, one concentrates 
on spectral functions Uj(Q 2 ): 

U 3 (Q 2 ) = i J d 4 xe^(0\T{j(x),j(0)}\0), (1) 
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where q 2 = — Q 2 and j(x) are local currents constructed on the quark and gluon fields. 
Then using the Operator Product Expansion (OPE) one obtains at large Q 2 : 

Tlj(Q 2 ) Hi Hj{Q 2 ) pa ,rton model ' (1 + ^ n Q2 Q*) ' ^ 

where the terms of order (InQ 2 ) -1 represent ordinary radiative corrections while bj are 
proportional to the (non-perturbative) gluon condensate, with the coefficient of propor- 
tionality depending on the current j. In principle, the fitted value of (a s (G^) 2 ) can 
depend on inclusion of higher orders in perturbation theory but the condensate appears 
to be large numerically and not sensitive to the subtractions. 

More recently, the gluon condensate was measured on the lattice. Here, one either 
subtracts the perturbation theory explicitly or keeps contribution of particular, pre- 
sumably dominating non-perturbative fluctuations ||. In what follows, we will always 
assume that the perturbative contribution to the condensates is subtracted. Moreover, 
we will propose a new method to obtain a lower bound on the condensates by considering 
their drop at the phase transition, see Sect. 5. 

A salient feature of (fj) is the absence of 1/Q 2 terms. An obvious candidate for 
d = 2 condensate is {(Affi). However, due to the gauge invariance the spectral functions 
Iij{Q 2 ) cannot depend on ((A®) 2 ) which is gauge dependent. The situation is changing 
if one considers a gauge non- invariant quantity, say, the gluon propagator Ii^ v {Q 2 ) j|, [|. 
Then the {(A a ^) 2 ) enters the OPE and its contribution in the £ gauge is: 

n^-d + fl|^^-f). (3) 

Since ({A a ^) 2 ) is gauge dependent it is not clear, however, what kind of information, if 
any, would be produced by measurements of ((A^) 2 ). 

For the framework which is outlined here, the crucial observation is that the minimal 
value of (A®) 2 may have physical meaning^ || [TJ. Indeed, consider a toy model when a 
plane is pierced by thin vortices carrying non-vanishing magnetic fluxes. Then, 

A»dx» = (flux) ^ 0, (4) 

where the integral is taken over a contour surrounding one of the vortices. It is clear then 
that (A 2 ) m i n is not zero and encodes information on the vortices. The idea that there 
exists a connection between ((AJ)^J and topological defects is the central one for the 
present review. However, the example of infinitely thin vortices in the continuum U(l) 
theory is in fact not adequate to represent topological defects since the corresponding 
action is infinite. By topological defects we will rather understand field configurations 
characterized by infinite potentials but surviving in the vacuum. 

In case of the compact U(l) theory the corresponding topological defects are Dirac 
strings with monopoles at the end points ||. The close connection between the ((Aj) min ) 

1 The content of the paper in Ref. M was summarized in the review article M . 
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and the topological defects was confirmed in this case through lattice simulations f| [/J, 
for a review see Sect. ||. 

Because of the asymptotic freedom, the topological defects in non-Abelian theories 
are associated with infinite potentials but with finite action. Since the topological defects 
are marked by potentials rather than by action, it is partly a matter of gauge fixing 
which topological defects are relevant to the vacuum state. Magnetic monopoles are 
apparently of special interest because of the dual-superconductor confinement mechanism. 
The magnetic monopoles do bring in a non-zero value of ((^)mm)- m Sect. |]we discuss 
the connection between the monopoles and ((^)mm) following the paper in Ref. |10| . 

An outcome of this consideration is a two-component picture for ((A^ in ), see Sect. [§ 
Namely, contributions of order Aq CD to ((A^ in ) are associated both with large (of order 
^qcd) an d small (of order a) distances, where a is the lattice spacing vanishing in the 
continuum limit. The soft part of ((^)min) can ^ e measured on the lattice by using 
the relations like (H). First measurements of this type were reported very recently jlT| 
(see also [|12j). By chance, the measurements of ((Affi) were performed in the Landau 
gauge, which is just the gauge which minimizes {A^) 2 . The measurements indicate that 
((A^)^ in ) is large numerically. We discuss phenomenological implications of this result in 
Sect. |5]. In this section we present also an independent estimate of ((Af l ) 2 nin ) by measuring 
its change at the deconfinement phase transition. 

It is most intriguing whether the "hard part" of ({A'J J ) 2 nin ) associated with the ultravi- 
olet region can enter any physical quantity. We speculate that the hard part of ((v4^)^ in ) 
could be related to the 1/Q 2 corrections discussed recently within various phenomenolog- 
ical frameworks (see [14], |l5j and references therein). Although such a relation cannot 
be proven, it could, in principle, be tested in the lattice simulations, see Sect. ||. 

2 Topological defects and ((^J^m) in the Abelian 
case. 

Let us consider first photodynamics, i.e. the theory with the Lagrangian 

L = ^(d,A u - d u A,) 2 . (5) 

Since there is no interaction at all, the coupling e 2 is not running. However, one can 
consider different theories with various values of e 2 . Although the Lagrangian looks trivial 
the physics depends actually on whether one considers non-compact or compact versions 
of (H|). It is only the non-compact formulation which describes free photons at any value 
of e 2 , while in the compact case there is a phase transition at e 2 ~ 1 which is due to the 
monopole condensation || (for a review see Hl6f). 

Physically, the difference between the compact and non-compact versions is that in 
the former case the Dirac strings carry no action and are allowed in the vacuum. The 
compact formulation naturally arises within the lattice regularization, where the action 
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is given in terms of the contour integrals around elementary plaquettes p: 



S = S p , S p = Re exp i (b A^dx^ . (6) 



( ■' p 



The standard continuum action @ emerges only in the limit of small gauge potentials 
after applying the Stokes theorem. However, the configurations for which 

A^dx" = 2nk , k e Z (7) 

for every plaquette evidently have vanishing action on the lattice. Thus the Dirac strings 
which are defined by Eq. (|7|) are the first example of what we call topological defects. 
These are field configurations which are constructed on the singular in the continuum 
limit gauge potentials but which do not cost infinite action. 

In fact, a closed Dirac string satisfying the condition (^) everywhere is not a proper 
example of topological defects. The point is that such a field configuration is to be 
considered as a gauge copy of the classically trivial vacuum, A^ = 0. In particular the 
closed Dirac strings would not contribute to ((^)mm) since it could be gauged back 
A^ = by means of the gauge transformations. Although the corresponding gauge 
transformations are singular in the continuum limit they should be admitted into the 



theory to match the lattice formulation, for more details see |[L7 |. 

The physically significant topological defects are therefore Dirac strings with open 
ends, or monopoles. Due to the magnetic flux conservation the action associated with the 
monopole is not vanishing at all but rather diverging in the ultraviolet: 

t r°° 

A ~ 4 / H 2 r 2 dr ~ ^ 



e 2 /„ e 2 a 



where L is the length of the monopole world-line and a is the lattice spacing which plays 
the role of the ultraviolet cut off. Because of (§) the monopoles are highly suppressed by 
the action for L ^> a. However, there is an entropy factor of the order e +( const ) L / a ^ where 
the constant (const) is of pure geometrical nature and is independent, of course, on e 2 . 
The entropy overcomes the suppression due to the action at e 2 ~ 1 and there is a phase 
transition, which is nothing else but the monopole condensation. 

Moreover, the model (|5|) provides the simplest example, where the dynamical relevance 
of ((^l/Jram) ma y b e shown analytically. In particular, the value of reflects the 

existence of topological defects (|7|||) in the theory. 

It is sufficient to consider one particular monopole trajectory. Therefore, we restrict 
ourselves to the non-compact version of (|5]). Then the monopole current j may be inserted 
into the vacuum via the 't Hooft loop operator []1 



H(A,Ej) = exp{S(F)-S(F + 27r*Z j )}, S(F) = j d 4 x(F^) 2 , (9) 
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where Ej is a Dirac string spanned on j and *E Mi , = ^^uXp^Xp- The insertion of the 
't Hooft loop is equivalent to the subtraction of the energy of one particular Dirac string 
Ej. Next consider the expectation value ((^4 M )min) m the theory 



Z{j) = J DA H(A, Ej) e~ s(F) = J VAex V {~^ J d 4 x(F^ + 2n*^ u ) 2 } . (10) 



Wmin) = - Jim - ^ lnZ y (i,m 2 )| m2=0 , (12) 



The gauge in which the A 2 is minimal is canonically fixed by the introducing the Faddeev- 
Popov unity: 

1 = A FP [\, A] J £>aexp{-A J d 4 x{A^ + d^af} (11) 

into the partition function (|l0|) and taking the limit A — > oo afterwards. Of course, this 
is equivalent to the direct fixation of the Landau gauge dA = 0. Then the expectation 
value ((AJ^J is to be calculated as: 

v-*x> V dm 2 

where Zv{j,m 2 ) is a partition function, defined in the finite volume V: 

Z v (j, m 2 ) = J VA 8{dA) exp{-^ J d 4 x [ {d {ll A v] + 2tt*E^) 2 + m 2 A\ } } . (13) 
The calculation of Zv{j,m 2 ) is straightforward and leads to: 

((AOLn) = Hm i ^[ const ■ indet(-<9 2 + 4e 2 m 2 ) + (14) 

+ V 2 jy -0 2 + 4e*m* + 47rW J v ^ S - -d 2 + 4e*m* ] ' 

where the first term, which is independent on the inserted monopole current, represents the 
perturbative contribution to ((A M )^ in ). Note that Eq. (|i~4D contains the string dependent 
term, which reflects the fact that A 2 is not a gauge invariant quantity. When calculated 
in a particular gauge A 2 depends on the position of the Dirac string. From Eq. (|14D we 
conclude that the non-perturbative part of ((.A M )j^ n ) 

C(e 2 ) = <(AJLn> - ((AJminW (15) 

depends only on the dynamics of monopoles and vanishes if these topological defects are 
absent. Thus ((e 2 ) should have a jump at the critical coupling although it cannot be an 
order parameter of the phase transition, since the monopole density is non zero even in 
the Coulomb phase. 

Fig. 1 represents the behavior of the quantity £(/3) as a function of the coupling 
constant (3 = 1/e 2 , calculated in the numerical simulations^ of the lattice compact U(l) 
gauge model. It clearly demonstrates that ((A^)^ n ) is indeed serves as a proper measure 
of the topological defects, at least in case of U(l). A unique and nice feature of the U(l) 
case is that the perturbative part of the condensate can be reliably removed. Indeed, 
in the non-compact case ({A^)^^) is entirely perturbative. Moreover, the perturbation 
theory is the same in the compact and non-compact cases. 

2 The details of the numerical calculation are not important for the present qualitative arguments. 
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Figure 1 (left): Nonperturbative part of ((A^)^^), Eq. ( |l5| ) in compact U{1). The phase 
transition occurs at = 1/e 2 as 1.0. 

Figure 2 (right): The densities of geometrical monopoles in SU{2) gluodynamics in the Landau 
gauge at zero and finite temperatures. The solid curves are the renormalization group prediction 



3 Geometrical monopoles and ((^)^ n ) 



Turn now to the gluodynamics, 



(16) 



where a = 1,2,3 is the color index (for simplicity we consider SU(2) gauge group only). 
If we would simply ignore the difference between the Abelian and non-Abelian cases 
and assume that the non-perturbative dynamics is the same, we would run into serious 
difficulties. Indeed, within the compact photodynamics ((^Omin) ~ a <1 since the UV cut 
off a is the only scale in the problem. In case of the gluodynamics there exists another 
(hidden) scale, Aq C d, defined in terms of the running coupling, g 2 (A.Q CD ) ~ 1. Moreover, 
naively the Abelian-like monopoles could not be relevant (non-perturbative) degrees of 
freedom in QCD since the action (||) is determined by the coupling at the UV cut off and 
g 2 (a~ 2 ) — > 0. As a result, the suppression of the Abelian monopoles due to their action 
is always stronger than the enhancement due to the entropy factor. 

The reality of the non-Abelian dynamics turns to be more varied. First, in the non- 
Abelian case the Dirac strings with open ends may cost no action at all [ I7| . The point is 



that although the Abelian magnetic flux is still conserved and the Abelian part of G a is 



singular at the end points of the Dirac string, the corresponding action is not necessarily 
large. Indeed, there is no direct connection any longer between the Abelian part of the 
field strength tensor and the action. The commutator term can cancel and does cancel 
the Abelian part in the explicit construction of Ref. [ PT| . 
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Thus, in gluodynamics there are no classical monopole-like solutions similar to the 
Abelian Dirac monopoles relevant to the U(l) case. This does not mean, however, that 
the monopoles are irrelevant in the non- Abelian case. To the contrary, the running of the 
coupling allows to scan the dynamics at various values of g 2 . Moreover, the same coupling 
governs dynamics of any U(l) subgroup of the SU(2). In particular, if g 2 ~ 1 then the 
condensation of Abelian-like monopoles is favored due to the entropy factor. And there 
is of course a lot of numerical evidence for the relevance of the monopoles, for review and 
further references see [fL{|. 

Here, we will outline the geometrical monopoles introduced in Ref. [fL0[| . Mathemati- 



cally, the basic idea is to define the U(l) subgroup relevant to the monopoles locally, for 
each plaquette. Physically, the crucial step is to relate the monopole properties to the 
minimization of (Afy 2 . 

We begin again with the observation that the action of SU(2) gluodynamics on the 
lattice is constructed in terms of elementary Wilson loops, not the continuum field strength 
tensor: 

U p = P exp i I A^dx" = e iFp = cos ^ + % n a p a a sin S , (17) 



F p = F p a a a /2, \F P \ = ^F«F*, (18) 

where a a are the Pauli matrices and we have defined n p = F p /\F p \. As a result, the field 
configurations for which 

\F P \ = 4irk, k e Z (19) 

cost no action and thus are the topological defects which we are looking for. Since the 
condition (|19|) is quite analogous to the Abelian case (cf. Eq. ([<])) we also refer to these 
topological defects as Dirac strings. In a closely related language, one can say that the 
Wilson loop (|T7D defines a "natural" U(l) associated with it as the group of rotations 
around the vector n p . In this way, one can define a U(l) group for each plaquette. The 
definition of the £7(1) subgroup varies from one plaquette to another, emphasizing the 
non- Abelian nature of the underlying theory. 

Knowledge of the Wilson loop (|T7|) does not allow, however, to determine k, Eq. ([19]). 
Indeed, the full plaquette matrix U p is the same for \F P \ = 4tt and \F P \ = 0. What is 
needed is the decomposition of the plaquette variable \F P \, 

\F p \ = ^1 + ^2 + ^3 + ^4, (20) 

in terms of the phases i = 1, ...,4, which are associated with the corresponding links. 
The decomposition ( p0|) comes about naturally in the basis of the coherent states (for 
details and further references see JL0|). Indeed, the Wilson loop W(T) is defined as an 
evolution operator of the quantum mechanical system, the state space of which carries 
the irreducible representation of the gauge group: 

(d t + iA)\1>) = 0, |V(T)> = W(T)\iP(0)). (21) 
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Eq. (|21|) implies that the evolving vector \ip(t)) is a generalized coherent state. Moreover, 
for a given evolution operator W(T) there always exists such a state | i/jq ) for which the 
entire evolution reduces to a phase factor: 

\MT)) = W(T)\M0)) = e^\M0))- (22) 

In terms of the phase <p(T) the fundamental Wilson loop is given by Tr W(T) = cos <p(T). 

In order to construct the decomposition (|20| ) we are making use of the following prop- 
erty of the coherent states: 

g\ij) = e^|V s >, VgeSU(2), (23) 



where | ifj g ) and (p g depend on both | ip ) and g. Then the decomposition fl2~0|) emerges as 
follows. For a given plaquette matrix U p one finds the eigenvector | i/; ): 

U p \^o) = e^l^o) (24) 

and then compute \F p \ using Eq. (p3|): 

C/ P | Vo > = f/i...f/ 4 |^o) = e*<*/*u 1 ...U 3 \fa) = e^ 1+ " + <^/ 2 | > (25) 

As a result, for any given lattice fields configuration, one can distinguish between \F p \ — 
Ank and \F p \ = for every plaquette and detect the Dirac strings in this way. The end 
points of the strings are then identified with monopoles. 

The construction outlined above fully determines the Dirac strings and monopoles as 
geometrical objects. As a mathematical construct, it certainly appears very appealing. 
However, from the physical point of view it is crucial that the monopoles constructed 
in this way are gauge dependent. Indeed, it is the decomposition (|20|) that determines 
whether a particular plaquette is pierced by a Dirac string. But the phases <pi are de- 
pendent on the link matrices and, therefore, gauge dependent. Moreover, contrary to the 
Abelian case even the end points of the strings (monopoles) are gauge dependent. 

We need at this point a physically motivated choice of the gauge. In view of our 
discussion, the Landau gauge which minimizes the (A^) 2 seems to be singled out. Indeed, 
in the continuum limit both the Dirac strings and monopoles correspond to singular gauge 
potentials. It is easy to imagine, therefore, that one can generate an arbitrary number of 
spurious singularities by going to arbitrary large potentials A a , so to say inflated by the 
gauge transformations. On the other hand, by minimizing {A^) 2 one may hope to squeeze 
the number of the topological defects to its minimum and these remaining objects may 
be physically significant. 

It is amusing that the guess on the choice of the gauge can be checked through nu- 
merical simulations. Indeed, if the geometrical monopoles are physical, then the lattice 
monopole density pi at should satisfy the renormalization group equation: 



Pphys 

plat = "4X3" 



153/121 , Qn2 



exp ( -—(3 ) , (26) 
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where p p h ys is the physical density, /3 = 4/g 2 and the scale parameter A is fixed by 
numerical value of string tension ay/a = 0.1326 at f3 = 2.6 (see, e.g., Ref. f2(|). The 



condition ( P6|) is a very strong constraint on the pi at and there is no surprise that if 
we do not fix the gauge in a particular way pi at does not satisfy (pop. However the 
geometrical monopoles defined in the Landau gauge turn to be physical objects, i.e. their 
density satisfies the condition ( p6| ) numerically. On the Fig. 2 we plot p\ a t in the Landau 
gauge versus (3 on the symmetric 12 4 and asymmetric 4 x 12 3 lattices (the latter case 
corresponds to a finite physical temperature). At zero temperature the lattice monopole 
density sharply follows Eq. fl26|) with \n(p phys / '4 A 3 ) rs 12.2, which is represented by the 
upper solid curve on the figure. Therefore, we can estimate the density of geometrical 
monopoles in 577(2) gluodynamics at zero temperature (confinement phase): 

P P t T) « (l-9v^) 3 « (840 MeV) 3 , (27) 

where the conventional value y/a = 440 MeV has been used. The monopole density in the 
high temperature deconfinement phase apparently scales in accord with Eq. fl26|), albeit 
with a somewhat smaller value of p p h ys '- 

Pjgt T) « (l-7v^) 3 « (760 MeV) 3 . (28) 

To summarize, the monopoles belong both to field theory and statistical physics. In- 
deed, monopoles are defined in field theoretical language. However, they are gauge depen- 
dent and any particular monopole can be removed by a gauge transformation. Therefore, 
the natural question is whether all of them can be removed by an appropriate choice of 
gauge. The gauge which maximally suppresses the monopoles is the Landau gauge de- 
fined by minimizing {A^) 2 . It turns out that precisely in this gauge the monopoles become 
physical. 



4 Anatomy of ((A^J 

Establishing the connection between the ((A°/)' 2 nin ) and the monopole physics allows to 
understand better the structure of the {{A^ l ) 2 nin ) itself. Indeed, both monopoles and Dirac 
strings contribute now to {{A a ^) 2 nin ). In this sense the {{A a ^) 2 nin ) condensate is basically 
different from, say, the gluon condensate. 

As for the monopole contribution, it is divergent in the infrared for a single monopole, 

j(AXond 4 x ~ L-R, (29) 

where L is the length of monopole world-line and R is the infrared cut off. In reality, it 
means that the contribution of the monopoles to {{A a ^) 2 nin ) is of order Aq CD . Indeed, at 
distances ~ ^qcd ^ ne approximation of the monopole gas is no longer valid. 
For a Dirac string we have an estimate: 

j{A a X trmg d'x ~ Hk uv ) ■ L ■ T , (30) 
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where huv ~ a ~ l is the ultraviolet cut off, L is the length and T is the time of existence 
of the string. The string is infinitely thin in the continuum limit and the contribution 
( |30D comes from singular potentials. There is a logarithmic divergence in the ultraviolet 
but this could be compensated if we consider (g 2 (Af l ) 2 nin ) instead of ((A° M )^ nin ). Then we 
can go to the continuum limit and, in the logarithmic approximation, there is no sign of 
the size of the string left. The contribution of the strings to ((A°)^ in ) is controlled by 
typical values of L, T. The typical values of L, T are not determined, however, by the 
dynamics of the strings themselves since the Dirac strings carry no action^. The Dirac 
strings emerge as a supplement to the monopoles, the dynamics of which is governed by 

9 2 (^qcd) ~ L That is wh Y L > T ~ A qcd as far as {( A lfmin) is concerned. 

It is worth emphasizing once more that no insight into the monopole dynamics can be 
gained through the quasiclassical approximation. In the classical approximation the Dirac 
string with open ends is not associated with any action |T7|]. One-loop corrections were 
also considered explicitly and do not distinguish this field configuration from the trivial 
vacuum either. Presumably, this is true to any finite order in the perturbation theory. 

So far we discussed connection of ((A*)^) with monopoles. However, the notion of 
((^)mm) may be more general than the monopole-related mechanism of confinement. 
Indeed, turn to another mechanism, that is P-vortices (for review and further references 
see, e.g., In this case the P-vortices seem also to be constructed on the topological 

defects, i.e. field configurations with singular potentials and no action at short distances. 
Thus, generically we would get again a two-component picture for {(A")^ nin ). Numerical 
studies along these lines would be very interesting. 



5 Measuring ((A^-J 



As is mentioned above, one can use relations like @ to determine {(A^) 2 ) from fits to 
the data. That is what was proposed some time ago f|, [5] and attempted very recently 



in the lattice simulations 11 



There are two comments on this approach which we would like to add. First, it is only 
the soft part of the ((A^) 2 ) which can be treated consistently via OPE. While in reality 
we expect that ((A a ^) 2 ) is contributed also by short distances, see above. Second, mea- 
surements in the Landau gauge are in fact singled out since then one measures not mere 
gauge artifacts but rather ((^)mm) which is physically meaningful. It is just happened 
so that the first measurements of {{A a ^) 2 ) have been performed in the Landau gauge and 
can be, therefore, interpreted in physical terms. 

In more detail, measurements of the 1/Q 2 corrections both to two- and three-point 
Green functions have been reported. In the latter case the measurements refer to the 
symmetrical point G(Q 2 ,Q 2 ,Q 2 ). The results of the fits [IT| are: 



(9 2 (AXin) = (2.32(6) GeV) 2 , (^K)LJ = (4.36(12) GeV ) 2 , (31) 



3 Here we discuss only monopoles and Dirac strings living in the vacuum. External monopoles can be 
introduced via the 't Hooft loop. The corresponding Dirac strings have in the continuum limit an infinite 
action, for further discussion see Jl7t|. 
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where the two numbers refer to the fits to 2- and 3-point Green functions, respectively. 
There is a discrepancy of factor about 4 between the two fits (|3l|) which might be 



due to the yet-inconsistent treatment of higher orders in perturbation theory [11]. Let us 



assume following Ref. [nj] that the account of the higher orders does not change the scale 
of the 1/Q 2 corrections. To appreciate the numbers (|3lf) it is convenient to introduce a 



tachyonic gluon mass WA 



(A a »(-q)Al( q )) = -^(<V - M^) « (1 - 5) • ?(<W - ^ ■ (32) 
p q z + g z q z q z q z 

Then the fits (|3lr) give m 2 ~ (1 -T- 4) GeV 2 . Moreover, since ((A®) 2 ) measured in the 
Landau gauge coincides with ((ijj^j which is physical, this m 2 g gives a measure of the 
non-perturbative corrections. And it is then for the first time so that numerically large 
1/Q 2 corrections have been directly observed. This result, in turn, is a confirmation of the 
theoretical speculations that there exists an "intermediate" mass scale which is formally 



of order Aq CD but is numerically large, see [§, 0, |2 

In view of the discrepancy between the two fits ( ]3"1"D it would be desirable to get 
an independent estimate of the ((^)mm)- m fact, a lower bound on ((^)^„) can 
be obtained in a remarkably simple way. Indeed, since the perturbative, or ultraviolet 
divergent part of ((^)?wm) * s n °t sensitive to the phase transition, the drop in ((A° i ) 2 nin ) 
at the critical temperature is due to a change in the non-perturbative part. Of course, 
the A 2 condensate might be non-zero in both phases, but since it cannot be negative 
the difference between ((i^ in ) above and below critical temperature provides a lower 
bound for {{Alf min ). 

We have considered pure SU(2) gauge theory on the 12 4 and 4 x 12 3 lattices in the 
Landau gauge and measured directly the quantity 



^^T(l-iTriyx)), (33) 



AV ^ 2 

x,ii 

where V denotes the lattice volume and U^x) are the link matrices. In the naive contin- 
uum limit (^) reduces to 

V - ^ " \ Jd'x ((AXJ (34) 

and includes both perturbative and non-perturbative contributions. The results of our 
measurements are summarized on the Fig. 3. It is straightforward then to estimate the 
drop in the function r](j3) across the deconfinement phase transition and obtain 



((AXin) > J^^V = = (761.6 MeV) 2 . (30) 



Numerically, the estimate (p5|) is in agreement with (|3T|). Moreover, Eqs. (|35|), (|27|), 
( p8|) allow us to speculate that the ((A'^) 2 nin ) condensate is mostly due to the monopoles. 
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Figure 3: The quantity 77, Eq. (33) versus (3 in SU(2) lattice gluodynamics. The solid curves 
are drawn to guide the eye. 



Indeed, the monopole density changes across the phase transition by approximately 20%, 
see Eqs. (pT|-[2"8|). If we naively scale the change in ((A°)^ n ) with the change in the 



monopole density we get estimates of the itself, which is quite close to (plQ . 



Since we treat the perturbative contributions very differently from the paper in Ref. ITT 
the result fl35|) can be considered as an independent confirmation of the high mass scale 
associated withg the ((A a ) 2 min ). 



6 Short-distance physics and ((^)^ n ) 

As is argued in Sect. 4, ((^)^ n ) is contributed also by short distances. This contribution 
cannot be treated by means of the OPE. The most intriguing question is whether this 
hard piece of ((A a ^) 2 min ) may have physical manifestations. 

The ((A" J ) 2 min ) / Q 2 corrections to the propagator can be conveniently traded for the 



gluon mass, see Eq. (|32|). In fact the notion of the gluon mass is more general than the 
Eq. d||) based on the OPE since the 1/Q 2 correction can be associated with short distances 
as well. Quite remarkably, the phenomenology of 1/Q 2 corrections to gauge invariant 



quantities in terms of a tachyonic gluon mass was proposed on heuristic grounds |L4|, 
and turned surprisingly successful. Moreover, there could be a close connection between 
the short- distance contribution to ((A't) 2 nin ) discussed above and the tachyonic mass. 
Indeed evaluating the A 2 condensate in terms of the propagator (ET) and subtracting the 



4 Note that we used the normalization of gauge potentials (see Eq. (jig)) different from that in Ref. ]TlJ 
and hence there is no explicit g 2 factor in (|3^). Moreover, the authors of Jll| considered SU(3) gauge 
group, not SU(2) as we did. However, the string tension is equal in both cases and, therefore, numerically 
the estimate (p3) should not differ considerably for SU(2) and SU(3). 
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perturbative part we get: 



<KDLn> = j^WHW ~ rn 2 g HA uv ) (36) 

and, therefore, m 2 conveniently parameterizes the sum of short- and large- distance con- 
tributions. Note that the tachyonic nature of the gluon mass reflects the positivity of 
the non-perturbative ((ip^ in ) in this approach while in Ref. (|TJJ the tachyonic sign of 
the m 2 was introduced on pure phenomenological grounds. In principle, one could hope 
that Eq. ( |3T)] ) can produce a kind of self-consistency equation. However, such an equation 
would be infrared sensitive and we we would not pursue this line of investigation here. 

Note that an infrared sensitive tachyonic gluon mass was introduced first by V.N. Gri- 
bov p3| . The physical meaning of this mass is that because of the hadronization gluons 
" decay" into hadrons. This is a pure non-perturbative effect and, within the perturbative 
expansion it is described as a "leakage" from the basis of the states used. On the theo- 
retical side, therefore, the central question is whether a similar interpretation is possible 
for the ultraviolet sensitive gluon mass which we are discussing now. The short- distance 
contribution to ((A°)^ in ) comes from field configurations which are given by singular po- 
tentials but possess no action. The Dirac strings are an example of such configurations. 
The question is, therefore, whether interaction of the gluons with the Dirac strings can 
be characterized by a tachyonic gluon mass. At first sight, it is hardly possible since the 
gluons should not interact with the Dirac strings. However, if one uses the perturbative 
basis of the plane wave functions the gluons do interact with the Dirac strings since the 
Dirac veto is not satisfied (for a recent review and further references see [pH]). In other 
words, accounting for the Dirac strings asks for a non-perturbative reshuffle of the basis 
and, within the perturbative approach, the effect might be described by a tachyonic gluon 
mass. But, of course, these considerations are highly speculative and we turn back to the 
phenomenology. 

Phenomenologically, the problem is whether one can separate, via measurements, the 
effect of ((i°)^ n ) entering through ordinary OPE from the effects of the short distance 
tachyonic gluon massQ. Physically, of course, the pictures are different since the effect 
of the soft part of ((A° x )^ nin ) cancels from the gauge invariant quantities while the effect 
of short distances may well persist. However, for a non-gauge invariant quantity the 
predictions can be the same or similar. In particular, Eq. (|32|) was postulated in f[4 



and m 2 was meant to parameterize the short- distance contribution to gauge invariant 
observables. Now, the propagator can be affected also by soft part of ((i°)^J, see 
Eq. (|3T)|). Moreover, Eqs. (j3"2D, (|3T)| ) are in fact identical. To distinguish between two 
contributions we will denote the short-distance part of m 2 as A 2 . Let us discuss the 
measurements of the 1/Q 2 corrections performed so far and their interpretation: 

5 In the quasiclassical approximation the non-perturbative power corrections from short distances are 
due to small-size instantons and suppressed by a high power of Q 2 , see Ref. |27| ]. This conclusion does 
not apply to our analysis of the monopole-related effects since monopoles cannot be obtained in the 
quasiclassical approximation, see above. 
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(i) To set the scale A 2 of the l/Q 2 corrections to gauge invariant quantities, let us men- 
tion that measurements of the heavy quark potential at short distances in gluodynamics 
indicate: 

A 2 ~ 1 GeV 2 . (37) 
This estimate obtained first in |TJ[ from the data on the full potential was later confirmed 



by analysis of the data on non-perturbative QQ potential at short distances [§] (see also 
|]28[| ). Note that in the realistic QCD, with the effect of the light quarks included the 
overall fit to the correlation functions (0) with inclusion of the l/Q 2 terms gives: 

A 2 ~ 0.5 GeV 2 . (38) 

(ii) The l/Q 2 corrections to the two- and three-point Green functions. As is already 
mentioned above, studying the corrections to the propagator does not allow to separate 
the short- and large- distance contributions to the gluon mass. The same is true in fact 



for the G(Q 2 ,Q 2 ,Q 2 ). As is noted in Ref. [O the effect of the soft part of ((^) 2 ) on 



the a s (Q 2 ) defined in terms of the 3-point Green function at the symmetrical point is 
determined by OPE: 

a s (Q 2 ) - aT rt (Q 2 ) ■ (1 - [f^l p ■ (39) 

Moreover the l/Q 2 correction in ( |39"D is entirely due to the renormalization of the external 
legs, or renormalization factor Z 3 ^ 2 . A similar equation holds for the corrections of order 
A 2 . Thus, introducing A 2 does not help to reduce the discrepancy between the two fits, 
see 

We pause here to note that it was not specified in Ref. WM in which gauge one intro- 



duces the short distance tachyonic gluon mass. As a result, the discussion was confined 
to one-gluon exchange. In view of the relation between ((^4 2 )) and the gluon mass, we 
would parameterize the l/Q 2 corrections in terms A 2 specifically in the Landau gauge. 

(m) The l/Q 2 corrections to the 3-point Green function at the asymmetrical point, 
G(Q 2 , Q 2 , 0). The corresponding measurements are reported in |L2| and the observed l/Q 2 
corrections are numerically large. The OPE does not apply in this case and, therefore, 
there are no predictions for the effect of the soft part of the ((A'p 2 nin ) [p!T| . As for the 
model where the whole effect is due to a short- distance gluon mass, the l/Q 2 corrections 
are due to the the same Z factors and the same as (|3"5|). (The persistence of the factor 
^•3/2 j g ^ ue £ p ar ^i cu i ar definitions of the three-point function accepted in [0). The 
data are indeed well fitted by A 2 ~ 1 GeV 2 . 

To summarize our discussions, all the l/Q 2 corrections observed so far could be ex- 
plained by the model 0, 24| with a tachyonic gluon mass A 2 ~ 1 GeV 2 . However, 
the quality of the data and their analysis is such that it is not ruled out at all that the 
gauge- dependent quantities, like the gluon propagator, receive also comparable contribu- 
tions from the soft part of ((A*)^ in ) which cancels from the OPE for gauge invariant 
quantities. Further measurements would hopefully clarify the situation. In particular, 
measurements of the 3-point function for all external momenta large but not equal would 
especially helpful. 
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Conclusions 



The minimal value of the potential squared, ((Ap^ jn ) encodes information on the topo- 
logical defects in gauge theories. Already first measurements of ^ n ^ ne com P ac t 
U{1) || 0] indicate that ((A'f l )^ nin ) is sensitive to field configurations responsible for the 
confinement. In that case these are monopoles 0, for a review see Sect. 2. 

Within the dual-superconductor mechanism, monopoles play also central role to ex- 
plain the confinement in QCD. Close connection between the ({A'J J )f nin ) and monopoles 
was revealed first in Ref. fT(| in terms of the so called geometrical monopoles, see Sect. [| 
Here, we extended the analysis of the numerical data on the density of the geometrical 
monopoles at temperatures below and above the deconfinement phase transition. Note 
that in the non-Abelian case it is partly a matter of gauge fixation, which non-perturbative 
fluctuations dominate in the infrared region. In particular, analysis of possible connection 
between the P- vortices and ((A")^ n ) is still awaiting its time. 

A novel feature of the ((i^ ffl ) is that even non-perturbatively it is contributed not 
only by large but small distances as well. The both contributions to ((A")^ in ) are of order 
A.q CD . The physics is that the density of topological defects is decided by interactions at 
large distances ~ Aq^ d . If one focuses at short distances, then the topological defects 
are build up on singular potentials which cost no non-Abelian action, however. And these 
singular potentials bring in a finite, up to logs, value of ((A^ in ). 

There are various ways to measure ((^)^ in ). First, one can attempt direct measure- 
ments on the lattice. The main problem here is the subtraction of the trivial perturbative 
contribution. The subtraction is easy to make in case of the compact U(l) || |7j, see 
Sect. 0. In the non-Abelian case we were able to establish a lower bound on the non- 
perturbative ((^)mj„) which is the drop in ((^)mm) across the phase transition, see 
Sect. §. 

The soft part of ((-A°)^ in ) enters OPE for gauge-variant quantities like the gluon 
propagator j|, |^, [IT], [12|| . Recently, numerically large l/Q 2 corrections were found to 
two- and three-point Green functions in the Landau gauge [TT], [EJ. Interpreted in terms 
of the ((^)mm) the data agree with the lower bound on ((A° M )'^ nin ) discussed above. 
Moreover, the large value of the l/Q 2 corrections confirms theoretical speculations that 
the actual scale of the violation of the asymptotic freedom can be large numerically, see, 
e.g., |H, [14], [15], and references therein. 

Separation of the short- and large- distance contributions to ((i°)^ in ) and, more gen- 
erally, to the l/Q 2 corrections remain a challenge to theory. The main problem is to clarify 
how the short distance contribution to {{A < f/) 2 niv ) enters various physical quantities. The 
hard part of ((A a ) 2 min ) could well be the fundamental structure behind the tachyonic mass 



introduced phenomenologically in [|1J] . If one accepts this assumption, the short distance 
contribution is sufficient numerically to explain all the existing data on the l/Q 2 correc- 
tions, within existing uncertainties of the analysis. A comparable contribution of the soft 
part of the ((vlj)^,,) to gauge-variant quantities (in the Landau gauge) is not ruled out 
either. 

In short, we believe that already now one can conclude that the dimension d = 2 
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\a\2 

allows for a new insight into the physics of both large and short distances in QCD. 



condensate in gauge theories, {(A'p^ nin ), encodes important dynamical information and 
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